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In order to explore the relationship between future observations of the precession of 
orbiting gyroscopes (GP-B) and existing observations of post-Newtonian, relativistic effects in 
gravity, we have formulated a phenomenological gravitational equation of motion for matter 
which presupposes the very minimum possible concerning the gravitational interaction. 

The motional part of the gravitational equation of motion — proportional to (v/c) 2 times 
the Newtonian accelerations — is responsible for the interaction of gyroscopes with other bodies, 
but also plays a role in many other observables in the solar system which have been measured 
already. The set of all observations is a redundant — completely determining but to date 
consistent — set, and they predict the gyroscope precession rates with better than ten percent 
accuracy. A number of anomalous gyroscope precession terms are calculated within the general 
model, but their existence would be in conflict with other non-gyroscope observations. 

' - ' -r:_ If the gyroscope precession rates are found to differ from the predictions which can be 

uniquely made from other observations of relativistic effects in the solar system, a major crisis 
and opportunity!) would face gravitational theory: it would then appear that no motional, post- 
Newtonian gravitational equation of motion can be constructed which reconciles the total set of 
post-Newtonian observations in the solar system. 


(NASA— CR— 187736 ) THE GRAV ITQHAGNF T IC N91-2la49 

INTERACTION ANO ITS RELATIONSHIP TO QTHFR 
RELATIVISTIC GRAV I T A T IGNAL EFFECTS Final 

Report, 1 Dec. 1990 - 30 Apr. 19 91 Unci as 

(Northwest Analysis) 54 p CSCL 03G G3/46 0003628 



NASA 

K*l«n» **•<«*/■( i aui 


1. Rtp<Kl No 


Report Documentation Page 


2. Government Accession No. 


4. TWt and Subtitle 

The Gravitomagnetic Interaction And Its 
Relationship To Other Relativistic Gravitational 
Effects 


//]/— 


S&'Z $ 


3. Recipient** Catalog No. 


6. Report Dete 

April, 1991 


6 Performing Of genii et ion Cod* 


7. Author (*) 


Kenneth Nordtvedt 


6. Performing Orgeniredon Nemo and Addreea 

Northwest Analysis 
118 Sourdough Ridge 
Bozeman, MT 59715 


12. Sponsoring Agency Nam* and Addraaa 


(406) 586-3263 


National Aeronautics and Space Administration 
George Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


8. Performing Organisation Raport No. 


10. Work Unit No 


11. Contract or Grant No. 

H-11215D 


13. Type of Report and Period Covered 

Final 12/1/90 to 4/30/91 


14. Sponaoring Agency Coda 


16. Supplementary Note* 

Modified version of this report to be submitted to Physical Review 
for publication* 


16. Abstract 

To better understand the relationship between the expected precession rates 
of an orbiting gyroscope (GP-B) and other observable consequences in the solar 
system of relativistic, post-Newtonian gravity, we have developed a phenomenological 
model of post-Newtonian gravity which presupposes the very minimum possible 
concerning the nature and foundations of the gravitational interaction. Solar system 

observations primarily interplanetary' ranging fix all the parameters in 

the phenomenological model to various levels of precision. This permits prediction 
of gyroscope precession rates to better than 10 percent accuracy. A number of 
new precession terms are calculated which would exist if gravity were not a metric 
field phenomenon, but this would clash with other empirical observations of 
post-Newtonian effects in gravity. It is shown that gravitomagnetism the 

post-Newtonian gravitational corrections to the interactions between moving matter 
plays a ubiquitous role in determining a wide variety of gravitational effects, 
including the precession of orbiting gyroscopes. 


17. Kay Word* (Suggested by Author! •» 

Gravitomagnetism, GP-B , relativistic 
effects in gravity 


18. Distribution Statement 

unclassified - unlimited 


19. Security Class#, (of thl* report) 

20. Security OatsJf. (of this page) 

none 

none 



22. Price 












I. Introduction 


Even before Einstein’s general relativity theory was created/ 0 there was informed 
speculation that the inertial properties of matter were determined by matter’s cosmological 
distribution, itself. (2) Such points of view became known as "Mach’s Principle", but often such 
concepts were stated without a specific quantitative model or theory for proper and complete 
formulation. 

It was soon noticed that in general relativity theory some of these Machian ideas were 
manifest. These concepts soon acquired the descriptive interpretation — "the dragging of inertial 
frames" due to proximate moving matter. 

In general relativity (and other metric theories of gravity) gravity is a consequence of the 
metric field g^fr.t). At any space-time location a coordinate transformation can be found that 
eliminates all first space and time partial derivatives of the metric field; and a simple linear 
transformation then can set g^ equal to special relativity’s Minkowski metric; 

gMv(r,t) V + 0(x f ' - x£) 2 U 

The locality is now a freely falling, inertial frame with only tidal gravitational fields. The 
transformation which accomplishes (1.1) is non-linear; 

* = (X- - o + )6 - xS) (X* - x£) 1.2 

I"J!p are the Chnsstofel symbols formed from g^ and its first partial derivatives, evaluated at the 
space-time locality x£ 

g 1 ™ (gav.3 + gpv.a " Sajr.v) 1.3 

For the space-time indices p, a, P taking various types of values, quasi-physical descriptions can 
be made of the various terms which occur in the transformation (1.2) to the locally freely falling 
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inertial frame. 


(p, a, P all spatial coordinates) These terms indicate a non-linear "warping" of the spatial 
coordinates. 

(p spatial, a and 3 time coordinates) These terms indicate an accelerative dragging of the inertial 
frame; 

6x k - Vi It, 5t* k = x,y,z 

(p spatial, a spatial and P temporal) These terms indicate a dragging of the inertial frame going 
linear in time; 

5x k - r£ 5x' 5t 

The Chrisstofel symbols T£ can be viewed as a 3x3 spatial matrix with dimensions of frequency 
(1/t). The antisymmetric part of then plays the special role of interest to this study; 

r ^ = s ik+ E e Uk w i 1.4 

i 

with s, k representing the symmetric part of rf c and the frequency vector to, being a convenient 
way to express the three independent components of the antisymmetric part of T*. e IIk is the 
antisymmetric permutation tensor. 

Whereas s, k represents a time dependent rescaling of the inertial coordinates (both a 
stretching and warping, in general), the interpretation of the effect of the antisymmetric part of 
r£ is unambiguous and familiar - a pure rotation of the local inertial frame; 

Sr(t) = - (0 x (r - r 0 ) 5t 

with © giving the rotation axis and rate of rotation. In weak gravity and the first post-Newtonian 
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order approximation; 

r\> “ ■ 8ito ■ goti) 1.5a 

with the antisymmetric part being simply; 

Ifco ' = gko.l * glo.k 1.5b 

The three components g to of the complete metric field are commonly called the gravitomagnetic 
vector potential in analogy with electromagnetic theory; 
gko ■ h k = x,y,z 

then 

© = - V x h i.(S 

The source of the gravitomagnetic vector potential is moving matter; in general relativity 
(in a common gauge); 

h = (7/2)^v(l/2)^v-ff i .7 

for a source m, moving at velocity v,. The acceleration of a test particle in the presence of a 
gravitomagnetic vector potential is 

fla = dh/ dt - vx(Vxty 1 g 

with dli/dt being an "electric" type acceleration and v x (V x E) being the "magnetic" type 
acceleration; 

6a(mag) = 411^/r 3 v x (r x v^ 19 
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Just as two electric currents produce a magnetic force between themselves, two mass currents 
produce a gravitational force (1.9) whose unique signature is that it is a gravitational acceleration 
proportional to the velocity of both interacting masses. (1.9) can be considered a different way 
of viewing the dragging of inertial frames. This interaction does not require spinning mass — 
simply moving mass. 

Superimposing the acceleration ( 1 .9) for a spinning body mass source and a spinning body 
test object, it then produces a torque acting on the two gyroscopes which causes precession of 
the gyroscope spin axes: from the geometrical point of view this precession of interacting 
gyroscopes is simply the following of rotating inertial frames. 

But, in fact, a gravitomagnetic contribution to the post-Newtonian gravitational interaction 
plays a ubiquitous role in contributing to observable phenomena.* 3 * It is present in any 
configuration of matter where there is mutual motion of two interacting bodies. The spin-spin 
interaction of two gyroscopes is just a special case and configuration for measuring the nature 
of this part of the complete post-Newtonian structure of gravity. It is not necessary to use 
gyroscopes in order to measure gravitomagnetism. 

There are other motional corrections to the static, Newtonian interaction. Just as ( 1 .9) is 
an acceleration proportional to (v)(v,), there are accelerations proportional to v 2 or The entire 
package of motional corrections will determine the properties of the gravitational interaction 
under Lorentz transformations. Gravity need not be Lorentz invariant.* 4 * If the underlying field 
theory of gravity were to be based on two (or more) tensor fields, or a vector field as well as a 
tensor field, then the cosmos could produce preferred inertial frames in its gravitational 
interaction: there would be phenomena resulting from the gravitational interaction which could 
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depend on the velocity of the system relative to the cosmological preferred frames. The 
gravitomagnetic interactions are part of the whole package of motional corrections to the static 
Newtonian interaction, which as a package determines whether gravity has preferred inertial 
frames or not. 

This study takes as phenomenological a point of view as possible. We build the post- 
Newtonian gravitational interaction structure from empirical observations, and with a minimum 

of a-priori theoretical presuppositions. Several assumptions which we will not make a-priori are 
that: 

1. Gravity is metric field based; 

2. Gravity has any special properties under Loientz transformations; 

3. Gravity fulfills any conservation laws or is derivable from a post-Newtonian many 
body Lagrangian. 

What is assumed is simply that there are post-Newtonian equations of motion for; 

1. mass particles in the presence of other matter; 

2. photons in proximity to matter; 

3. non-gravitational clocks near matter. 

Isolated, but general N body gravitational systems are assumed to be surrounded by 
asymptotic inertial frames, in which non-gravitational rulers and clocks fulfill the transformation 
laws of special relativity. The coordinate speed of light is assumed to be a function of space- 
time position relative to matter, and therefore it can be expressed by a phenomenological 
expansion; 
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1.10 



The dynamical evolution of the light propagation vector c can be based on (1.10) and a "least 
time" principle, or it can be assumed to be based on an independent phenomenological expansion; 

de/ds = r 2 £ + 1.11 


Since there have been independent observations in the solar system of the effect of proximate 
matter on light propagation times (1.10) <5) and on deflection of light (1.1 1)« we can empirically 
conclude that; 

r 2 = - r, - 2 ± 10 - 3 1.12a 

and 


Experiments on the rate at which non-gravitational clocks run when in proximity to matter 
(and in motion) find first that different types of non-gravitational clocks all behave in a universal 
fashion; (T) and secondly that in lowest order (in powers of 1/c 2 ) their proper rate t is given by: (i) 


dx/dt = 1 - a £ m,/|r -r, | - Vfev 2 + ... 


1.13 


with 

a = 1±i0 r 4 ! 14 

and t is the time kept by a non-gravitational clock which is at rest in the inertial frame and 
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located asymptotically far from matter. 

Non-gravitational clocks and photons, governed by the empirically determined expressions 
(1.10), (1.1 1) and (1.13), are then the chief probes with which to look at the dynamics of bodies 
in the solar system. It remains to assume a most general, post-Newtonian gravitational equation 
of motion for an N body system of particles. The coordinate accelerations of particle i in the 
presence of other bodies jjc will be assumed to be of the form; 

d2 ?/dt 2 = £ Gm^r 3 + £ gy(motional) + £ g^non -linear) (1.15) 

J i jk 

The first term in (1.15) is the Newtonian approximation to the gravitational interaction; the two 
post-Newtonian contributions to (1.15) are proportional to 1/c 2 : they differ in that the two body 
motional terms are proportional to (v/c) J — the square of body velocities — times Newtonian 
accelerations, while the non-linear accelerative terms in (1.15) have an extra power of Newtonian 
gravitational potential (Gm/c 2 r) multiplying Newtonian acceleration. 

The dynamical equation of motion (1.15) can only be assumed to hold in one 
cosmological inertial frame. Until the detailed structure of ft/motional) is specified by empirical 
observation, it can not be assumed that (1.15) has any special Lorentz transformation properties. 

As seen in (2.1) ^(motional) has ten arbitrary coefficients in its general expression, of 
which one is a free gauge or coordinate system parameter. This leaves nine coefficients to be 
specified by nine (or more) independent empirical observations concerning the shape and rate of 
solar system body orbits. If more than nine independent observations are available, but which 
depend only on the nine free motional coefficients, redundant constraint will exist on the model 
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(1.15). These redundant observations could conceivably clash and not all be consistent with any 
choice of the nine coefficients in the expression for ^(motional). We would view this as a 
fundamental crisis (and perhaps opportunity!) in gravitational theory, as it is difficult to see how 

any theory could be constructed which does not lead to motional post-Newtonian gravitational 
accelerations of the form (2.1) 

In particular, consider the future plans to orbit a precision gyroscope and measure the 
precession of its spin axis with respect to a reference frame defined by lines of sight to the 

distant stars. <9) There are expected secular precessions of the gyroscope’s axis due to two 
sources: 

1. It s orbital motion about the Earth (geodetic precession); 

2. It’s interaction with the gravitomagnetic field of the spinning Earth (inertial frame 
dragging). 

From the point of view of our phenomenological gravitational equation of motion (1.15), both 
types of gyroscopic precession are simply additional manifestations of the motional, post- 
Newtonian corrections to Newtonian gravity: the geodetic precession is viewed as a consequence 
of acceleration terms proportional to the square of the velocity of the test body mass elements, 
while the precession related to the dragging of inertial frames is considered a result of 
acceleration terms proportional to the velocity of source masses and the velocity of the test body 
masses. We will derive expressions for the general secular gyroscope precessions which can 
occur in an equation of motion of the type (1.15) (see actually (2.1)); and also we will obtain the 
many other observables which result from ^(motional) but are not gyroscope observables. There 
will be more independent observables than there are free coefficients in (1.15) — redundancy! - 
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- and therefore the possibility of inconsistency of gyroscope observations with the other solar 
system observations. 

There are also some other non-secular (but nevertheless possibly measurable) gyroscope 
precessions which result from the phenomenological model (1.15). Some of these non-secular 
precessions would imply non-metric gravity. The solar system is traveling relative to the cosmos 
at a speed of order |w| - 2 10 3 c, while an earth orbiting gyroscope is carried around the Sun 

at speed |v| - 10“*c. The equation of motion (2.1) with its w-dependent terms (2.3) can 

therefore possibly produce additional perturbations of the spin axis of a gyroscope which are 
dimensionally of the form; 

Cs = gjS-vv 1.16a 

and 

8s = gj •wv +g 3 s -vw +g 4 sw -v 1.16b 

and in which the coefficients g, , ... g 4 would be linear combinations of the phenomenological 
coefficients c, in (2.1). The term (1.16a) is only a few milliarcseconds in magnitude; however 
the preferred frame terms (1.16b) vary about 80 milliarcseconds with an annual period. 

The format of this study is as follows. In section II. we state the general 
phenomenological form that the motional, post-Newtonian gravitational acceleration terms can 
have in an N body gravitational system. Then a set of "observables" are derived which are 
attributes of orbits and dynamics of planets and test bodies in the solar system which are 
routinely measured by ranging experiments and/or telescopic observations. The observables form 
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a redundant set — they more than uniquely determine the free coefficients in the equation of 
motion. The observables could, in fact, possibly be inconsistent with the general model. 
Precession rates for orbiting gyroscopes are calculated as some of the observables. All possible 
gyroscope precessions are determined by coefficients measured in other non-gyroscope 
observations. 

In section HI. additional observations are derived which depend also on the non-linear 
structure of the post-Newtonian gravitational interaction. Some of these observations produce 
additional constraint on the motional coefficients c, ... c l0 . 

In section IV. we review the accuracy with which the various observables have been 
measured by past solar system observations, and how well this then permits us to predict the 
magnitude of the gyroscope precession terms. 

Section V. discusses some of the conservation laws in physical theory. The existence of 
conservation laws is shown to be related to constraints on the structure of the gravitational 
equation of motion, i.e. it requires relationships to hold among the coefficients c, in these 
equations, and therefore is predictive concerning some of the observables of gravity. 

Section VI. derives the constraints on the gravitational equation of motion which indicate 
that post-Newtonian gravity is derivable from a metric field. If gravity is a metric field based 
interaction, then several of the possible gyroscope precession terms vanish as a consequence of 
this principle, alone. 

In section VII. it is shown that there is a family of coordinate transformations which 
makes one of the ten coefficients in the gravitational equation of motion arbitrary. But it is 
shown how all ten coefficients change under this transformation, and this permits confirmation 
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that all observables calculated in this study are invariant under this coordinate transformation. 

Finally, in section Vin. we collect all the gyroscope precession terms derived in our 
general model of gravity. We discuss their magnitude, their relationship to the structure of the 
gravitational interaction, and the implications for theory if the observed gyroscope precessions 
were to clash with the predictions. 

II. Motional Post-Newtonian Corrections to the Gravitational Interaction. 

Since our goal is to specify the post-Newtonian gravitational interaction solely from 
experimental observations - including null observations - the starting point for analysis of 
observation is a general, phenomenological expression for the equation of motion for N particles 

(bodies of negligible gravitational self-energy or binding energy) which gravitationally interact 
with each other. 

A-priori this equation of motion can be assumed to hold in only one special asymptotic 
(r — ► o®) inertial frame (the preferred frame); Lorentz invariance of the gravitational interaction 
is not presumed but rather is to be discovered, if it exists, from the experimental observations. 
A system of bodies which moves as a whole when viewed in the special inertial frame will be 
analyzed by the N-body equation of motion, and the dynamical consequences of the system’s 
motion determined for comparison with observation. The solar system, for example, moves 
relative to the universe’s local preferred inertial frame at a speed of order 10' 3 c, yet no orbital 
effects related to this motion are evident within experimental accuracy. 

Post-Newtonian corrections to the Newtonian acceleration are of order 1/c 1 and take two 
forms. There are motional corrections which are of strength (v / cf relative to the Newtonian 
acceleration. The second type of post-Newtonian corrections are the nnn-li™»«r terms 


11 


proportional to the square of mass sources ( 0 ^ 01 ^) and therefore dimensionally proportional to the 
inverse cube power of interbody lengths (1/r 1 ). 

Writing a most general expression for the motional corrections, the phenomenological 
post-Newtonian equation of motion takes the form; 





tv 


V. + C S V , V J + C 9 V J V I + 


lj 


c .oVj] + £ m , m k Pijk 
jk 


2.1 


Units of G = c = 1 are used, p,* is a general vector expression of dimensions 1/r 3 and is 
composed from the inner-body vectors t i} , r Jk and r u . Consideration of the non-linear part of the 
post-Newtonain interaction will take place in section m of this paper. c„ c, ... c I0 are ten 
dimensionless coefficients expressing the most general motional structure of the post-Newtonian 
interaction. The coordinate (or variable) freedom of the equation of motion discussed in section 
VII indicates that a "gauge" can be chosen to eliminate any one of the ten c,. Our choice will 
be to set c 4 = 0; however this calculation of observables in this section will include all ten c,; this 
explicitly reveals the gauge or coordinate invariance of the calculated observables. 

In this section, experimental observables dependent only on ten c, of the motional 
structure of the equations of motion will be derived; section m will derive additional observables 
also dependent on the non-linear structure of the equation of motion. 

A system of bodies collectively moving at velocity w relative to the preferred inertial 
frame in which (2.1) is valid is considered. The coordinate position of each body then takes the 
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form; 


r, — > r, + wt 


2.2a 


and 


r i * v i -» v, + w (2.2b) 

For such a system the equations of motion given by (2.1) develop additional terms proportional 
to w and (w) 2 ; 


®I S "E —Ip f( 2c , + c 2 )w *v ( + (2c, + c 2 )w -V 

1 r <i 

+ (2c 4 + c 3 )w • fyV, • i tj - (2c 6 * c,)w -f (j ] 

_ m _ 

+ L, —r ’ [ (c 7 + W" + ( c 7 + c^wv, + (C, + c l0 ) Vj w + (c, + c 10 )™J 2.3 

J r »j 


-r ^[<c, 


+ C, 


+ c 3 )w 2 - (c 4 + c 5 + c,)(w -tyf] + (c 7 


+ C . +C 9 + C I0 


>r 

J 


m jV 


ww 


Post-Newtonian orbital effects in a system moving in the preferred inertial frame therefore 
naturally divide into w and (w) 2 effects; the latter are considered first. 

(Observable 1) Consider a two body orbiting system which moves in a direction w perpendicular 
to the orbital plane defined by r,j and v (j ■ dr tJ /dt. 

■i " I 1 + (c > + c 2 + c 3> w2 ] 2.4 

r !2 
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and a similar equation for body (2). The solutions to (2.4) are simply the Keplerian orbits but 
with periods altered by the factor; 


T(w) = T(o) 


1 ' i (c ' * c > 


c,)w 2 


2.5 


Atomic clocks far from the system and moving along with the system (i.e., clocks in system’s 
rest frame) will be running slow by the special relativistic factor (1 - Vi w 2 ). If gravitational 
clocks (orbits) are to show no dependence of their rate on their motion relative to the cosmos, 
the following constraint must hold: 


c, + c 2 + c 3 = -1 2.6 

(Observable 2) Now let the velocity vector w lie in the orbital plane. The equation of motion 
is now 


m,r , r , 

— [1 + (c, + c 2 + c 3 )w 2 + (c 4 + c 3 + c 6 )(w • t l2 f\ 

r« 

*ww 


2.7 


+ (c 7 + c, + c 9 + c 10 ). 


1 12 


Applying this to a nominal circular orbit or radius r c , there results a radial perturbation 


r w 


5r(t) = --y- [c 7 + c, + c 9 + c l0 - c 4 + c 3 + c s )]cos2tot 


2.8a 


But an observer at rest in die cosmos must "see" or measure a Lorentz contraction of the orbit 
Sfrt)^. = -I4r o w 2 cos2cot 2.8b 

Otherwise a co-moving distant observer will see contractions of the orbit proportional to the 
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square of the speed of the orbit through the cosmos. Comparing (2.8a,b) leads to the constraint; 

c 7 + c, + c, + c 10 - Vi(c 4 + c 5 + c 6 ) = 3/4 2.9 

(Observable 3) The period of this preceding orbit can be calculated; 


r 2 n "|2 rn, + nij 

1 + (c, + Cj + c 3 )w J 

|T(w>J 

+ 1/2(c 4 + c 3 + c 3 - c 7 - c, - c 9 - c 10 )w l 


But r 0 is the mean radius of the Lorentz contracted orbit; the transverse radius of the orbit is 
(using (2.8a)); 


, , 2(c 7 +c,+ c 9 +c, 0 )-c 4 -c 3 -c 6 2 


2.11 


Rewriting (2.10) in terms of r x yields; 


r T 

m, + nij 

, 2 ( C 1 + C 2 + C 3 ) + c 7 + c g + c 9 + c 10 2 

|_T(w) J 

3 

Tx 

2 


For these gravitational clocks to be properly time dilated as well, a constraint exists; 

C, + C 2 + c 3 + V4(c 7 + c 8 + c, + c 10 ) = -1 2.13 

A very accurate way to observationally test for the independence of orbital period dilation on the 
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orientation of the orbit relative to the cosmos is to measure the isotropy of the Newtonian 
interaction in celestial bodies moving in the cosmos. This will be discussed further in the section 
on empirical observations. From (2.6) and (2.13) the condition for isotropy of the strength of 
the Newtonian interaction is just; 

C7 + + c, + c 10 = 0 2.14 

(Observable 4) Consider the time rate of change of the angular motion for a two body system; 


d,- - m, - m, 

— ( r i2 x v , 2 ) = (c 7 + c 9 ) r 12 • w r I2 x v, 2 


1 12 


(c 7 + c 8 ) 


- - - 


r , 2 - v i 2 r, 2 xw 


‘12 


2.15 


The first term is a periodic perturbation, but the second term gives a secular perturbation to 

angular motion for eccentric orbits. The failure to observe such orbital effects requires the 
constraint: 

c 7 + c 8 = 0 2.I6 

(Observable 5) Consider the acceleration of a two body system’s center of mass; 


d 2 

_ _ \ 
"Vi + m2 r 2 

m ni _ 
^ - 1 2 r 

2(c, - c 3 )w • v 

dt 2 

m +m_ 
v 1 ^ J 

mr 3 

+ 2(c 4 - c 6 )w • fv • f 


+ m,m 2 (c 7 + c 8 -c, -c 10 jr -vw 
mr 3 + (c 7 +c , -c, -c |0 )r -wv 


For a circular orbit the projection of the center of mass acceleration in the direction w x (r x v) 
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is; 


• ^ a \ d 2 R m.rrLwv 

wx(fxO)' = - 1 2 

dt 2 


mr 


2(c, - c 3 )(w *9) 2 

|_ + < C 7 +c 9 -C* "C |0 )(W f) 2 


2.18 


Both terms in (2.18) produce secular accelerations; they cancel however if there is the constraint; 

2(c, - c 3 ) + Cj + c 9 - c g - c 10 = 0 2.19 

This constraint results in elimination of secular center of mass acceleration for eccentric orbits 
as well. The other two terms in (2.17) produce periodic accelerations of the center of mass. 
(Observable 6) The complete equation of motion for the interbody vector 7 = 7, - 7 2 of a two 
body system is obtained from (2.3); including the w terms it is; 



mr _ - m, 


(2c, + c 2 )w • v 
+ (2c 4 + c 3 )w -ff -v 


2.20 


ntj - m, _ 

+ p — [( c 7 + c *> r • VW + (c 7 + c,)r • wv] 


The non-secular contribution to angular motion variation is: 

d /- - , m, - m, 

— ( r x v) = (c 7 + c 9 ) L r-wrxv 

dt 


The equation for radial perturbation x(t) of a circular orbit is then; 
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X + co 2 x 


m, - m, 

— - — w *v [2c, + c 2 + 2 c 7 + 2c 9 jsinctX 


2.22 


t = 0 corresponds to r(t) traversing the in-plane component of w (w* is the magnitude of w lying 
in the orbital plane). (2.22) has a runaway solution — x(t) - tcosox ~ contrary to observation, 
requiring the constraint; 

2c, + c 2 +2(c 7 + c,) = 0 2 23 

(Observable 7) The angular motion equation (2.21) then results in an angular position 
perturbation; 

m, - nij 

ro0(t) = r(c 7 + c 9 )wv cosox 2.24 


For this longitudinal perturbation to be unseen after a Lorentz transformation to the inertial frame 
of a co-moving observer, a constraint is required; 


c 7 + c 9 = 1 

(Observable 8) From (2.20) the total radial perturbation is; 


2.25 


5a t 



(2c, + c 2 )w *v 

+ (2c 4 + c 3 - 2c 7 - c g - c 9 )w • ft • v 


2.26 


This radial perturbation along with the angular perturbation given by (2.21) can be applied to the 
general Keplerian elliptical orbit. The orbital perturbation is obtained most easily in the u(0) = 
l/r(0) representation of the orbit; 
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w(2c 4 + Cj)sin(20 - v) 


2.27 


5u _ e ( m i " n^) 
u o ~6 7 


with the Keplerian orbit being given by; 

1/r = u = u 0 { 1 + e cos0} 
and 

f = F x vl 

The angle v is the orbital longitude of the in-plane component of w. Failure to see such orbital 
bulges requires the constraint; 

2c 4 + c, = 0 2.28 

(Observable 9) Consider a three body system consisting of a massive central body of mass M, 
with another body of mass m in circular orbit. A third body orbits the body of mass m. From 
(2.1) there will be a difference between the acceleration of the test body m by the massive central 
body, because of the difference in the two bodies’ velocities. Collecting the terms linear in u, 
the velocity of the test body relative to body m, the relative acceleration is; 

5a = -2c -^Lv‘u+ c.-^Lr-uV 2 29 

1 R 3 7 r3 zzy 


R is the position of body m from the central body M, V = dR/dt and u= dr/dt; r is the position 
of the test body from body m. The secular part of the acceleration which results from (2.29) is; 
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Si = - 


\ 


— + c, 


J 


M- T - 
— u x L 
R 2 


2.30 


with L being the angular motion of body m; 


E= Rx V 

Applying this perturbation to a near-circular test body orbit which lies in the plane of the circular 

orbit of body m, the orbital frequency © and natural frequency for radial oscillations (eccentric 
motion) ©„ are altered; 




\ 


+ c. 


) 


M .. 
— V© 

R 2 


2.31a 


©„ 





2.31b 


The difference produced between these frequencies is a contribution to the precession rate of the 


test body’s orbital periastron; 


5(© o - ©) = -Vi (c, + Vi c 7 ) MV/R 2 


2.32 
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This is commonly called the "geodetic precession" of an orbit. Observations of the geodetic 
precession contribution to the Moon’s orbit produced by the motion of the Earth-Moon system 
about the sun produces a weak constraint; 00 ' 

c, + Vic 7 = 3 2.33 

(Observable 10) A test body is in a circular orbit about a body of mass m which gravitationally 
free falls toward another body M (M » m) which is at rest; the test body’s orbit around m is 
perpendicular to the direction toward M (R • 7 = 0). The equation of motion for the test body 
can then be written as; 


a 



1 +(c, 


+ c 2 +c 3 )V 2 ] + c 7 ^lR-ViI 

R 


MR + M mr 
~R r * X ~R~ 


2.34 


with V = dr/dt, u = dr/dt: R is the position of m relative to M, and r is the position of the test 
body relative to m. The x term is a necessary non-linear acceleration term which will drop out 

of the final expression for an observable. The gravitational fteefall of m fulfills the Newtonian 
energy constraint; 

ViV 1 - M/R = const 

The test body’s angular motion can then be obtained from (2.34); it evolves in time during the 
freefall; 

|fxu| . I - (_|l 2.35 


This, in turn, produces an evolution in the radius and angular frequency of the test body’s orbit; 
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r « r 

M 


2.36a 


1 + (X “ 2 c 7 - 2(c, + c 2 




= 6 


1 + ( 3c t ~ 2 X + 4 (c, 



2.36b 


in which the relationship V 2 = 2M/R was used. Atomic clocks freefalling in the gravitational 
field of body M will change their ticking rate according to; (,) 


© = © 


1 -(1 + a) 


M 


2.37 


while the proper radius of the orbit is given by 


r 

P 


1 - (a + D— 
R 


J 


2.38 


In (2.37) and (2.38) we have introduced two additional phenomenological parameters related to 
non-gravitational clock behavior and the light propagation rate in a gravitational environment; 

± = 1 -V*v 2 -aU(?) 2.39 


and 


c(r) = 1 + rU(r) 


2.40 


Observationally a = 1 ± 10 “* and T = -2 ± 10' 3 . 

A combination of (2.36a, b) can be formed to produce an observable dependent only on 
the motional coefficient c 7 along with coefficients related to the gravitational influences on light 
propagation and non-gravitational clocks; 


f \ 1 

f \ 

r 



—L - 1 

+ 

— L - 1 

r 

r~ J 


K J 




2.41 


in which the proper frequencies % are measured by co-falling nongravitational clocks, and the 
proper orbital radius is measured in terms of locally measured light propagation times across the 
orbit. In the solar system no variation of satellite orbital radii or periods are observed as these 
satellites follow their planets on the eccentric motion about the Sun: this requires the constraint: 
c 7 = 1 - oc - 2r 2 42 

A constraint on non-linear aspects of the equation of motion follows from the separate absences 
of anomalous size or frequency evolution of the orbits; 


( \ 


( \ 

r 


6 

— - 1 

s 

— - - 1 

r 

V ^ ) 


j 


= 0 


2.43 


requires the constraint; 
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x = 2(C, + c 2 + Cj) + 


2.44 


- C 7 + 

2 2 

But this constraint will be considered along with the other constraints on non-linear aspects of 
the gravitational interaction in section in (x = c,, + c (2 in the notation of that section). 
(Observable 11) Consider an extended body in motion at velocity w. The matter distribution in 
the body will be Lorentz contracted with respect to the proper coordinates, the contraction being; 


5rj = -l^Tj • ww 

This alters the Newtonian potential of that body; 

And the gradient of (2.46) in an altered acceleration field; 
8a * V8U = ViVw*Vw'^p J 7 


2.45 


2.46 


From (2.3) there is one term which generates a similar type acceleration field; 
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2.48 


5a = -(c 4+ c 3+ c 6 )£ — 4-(w-(7-r)) 2 (7-r) 
j |r - r t 5 J/ ' 

L ..... - 

-(e 4 * c 3 + c 6 ) — (w • i) 2 r - Vi Vw • Vw • ^ j. ^ r t 


For a spherical body 


y _ m *_ ? = I_L 

^ |r-r.| J 3 r 3 


2.49 


with 


1 " E m /j 2 


The anomalous acceleration field from (2.47) plus (2.48) is then; 


51 = ifi + c + c + c )± j l5 ( w ~ 6w twI 

9 2 4 5 6 r 3 -3W 2 ? 


2.50 


Absence of observed perturbations of low satellite orbits or gravimeter anomalies on the Earth 
surface requires the constraint; 

c 4 + Cj + c^ = -3/2 2.51 

(Observable 12) Let the extended body be spinning; 


V, = Qxr, 


2.52 
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Because of the time component of the Lorentz transformation, at a simultaneous preferred frame 
time there is unequal proper time across the extended body; 


Tj = t - W • Tj 


The rotational motion (2.52) along with the time transformation (2.53) then 
displacement of each mass element relative to its proper position; 

Sfj = V/Tj - t) = -VjW ■l j 

which alters the Newtonian potential of the body by; 


m, 


5U -- V ^7T^ 5F . 


nv - 


-v. E ^Tn^.7, 

and the acceleration field by; 


5a = V5U 


A number of terms in (2.3) produce a similar type of acceleration field; 

<•— m.(r — r.) _ _ 

5a = E |7 _ 7| , [< 2c j + c > • v j + (c 5 + 2c> • (t - tp Vj • a - y] 
+ E F 'V * + c .o) + <r - rp -wv/c, + C (0 )] 


2.53 
produces a 

2.54 

2.55 


2.56 
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with 


^ ra i Qxr i .«Jx7 


T being the body’s rotational angular momentum. The total anomalous acceleration field from 
(2.55) and (2.56) is then; 


5a = Vi 

( 


( , V 

c, + 2c. 

1 + c + 2c + 5 


Vi 


c « +c .o- 


C J + 2C « 


3 

V 


J 


wxj _ ^ 7 • w x j7 


V 


w x J r xJr • w 


2.57 


Absence of low satellite orbit perturbations or gravimeter anomalies on Earth which would be 
produced by (2.57) require the constraints; 

1 + c 2 + 2c 3 + (c 3 + 2c s )/3 = 0 2.58a 

c* + c 10 - (c 5 + 2c 6 )/3 = 0 2.58b 

(Observable 13) Consider the interaction terms in (2.1) proportional to (v,Vj) and the source 
consisting of a rotating extended body; 

v * O x 7j 

From (2.1) this source produces the acceleration field; 
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2.59 


C - C„ V 


a i = 


i xJ 3 ( 


2 


c. “ — 1 c, 

* 3 


_ V 

T.xjT.-v, 


c + c 

2 3 5 


i Vi x J • r , r i 


r, is the position of body i from the center of the source, J is the total rotational angular 
momentum of the source 


J = 53 m ]fj x (^ x V 

j 


2.60 


For an orbiting test body the time rate of change of the angular motion of the orbit is; 


dt 1 1 4 r J 


1 


C 2- C S + 3- C 5 


2 3r,- J v 1+ v t - J 7 t 


2.61 


which results in a secular precession of the orbit about T at the rate; 


(0 


C 2 +C . 


a 0 (l - e z ) 


2\3/2 


2.62 


a„ and e being the orbit’s semimajor axis and eccentricity, respectively. Such a precession is one 
manifestation of the so-called "dragging of inertial frames" by a rotating mass with angular 
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momentum | 


(Observable 14) The acceleration field (2.59) produced by one spinning mass causes a second 
spinning mass to precess. Let the test "gyroscope" have its mass elements moving with 
velocities; 


with the body having total rotational angular momentum; 
S * £ m t r ( x(f2xr,) 

i 


2.64 


The precession of 5 is then governed by; 


S 



c 2 +c « JxS J-RRxS 

8 f R 5- R s _ 

„ A 2 'j 7 x RR • S + J x R • $R 
c. + _ c, 

3 3 J R 3 


2.65 


The first term in (2.65) gives an instantaneous precession rate; 

- . c 2 + c, j-3j# 

8 R 3 


2.66 


while the second term produces a secular contribution to 5 of: 


S 




Jxfyp-s+Jx^-sp 

a„ 3 (l - e J ) J * 


2.67a 
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A 

p is the orbit s normal vector. (2.67) can have a non-zero component along 5 which would 
' change the magnitude of S. There are also a few non-secular (annual period) contributions to a 
gyroscope’s precession which are proportional to the velocity w of the solar system with respect 
to the preferred inertial frame. Measured in the solar system rest frame the variation of a 
gyroscope’s angular motion is given by: 


5S = Vz(4 - 3c, - c g - 2c 9 + 2cj + c 2 + c 4 + c,) V • w S 

S Vw 


7 ■ v g — 2 + 2c. + c 7 - 

1 # C 1 2 c .tf ... 


C 7 +C 

+ — -S wV 


2.67b 


plus the additional precession rate which is not integrable in closed form; 


dS/dt = — — Rwft&S 

2 R 3 


2.67c 


These precessions have dimensional amplitude of 40 milliarcseconds (mas) for motion around 
the Sun on the Earth’s orbit. But all these expressions (2.67b,c) have zero coefficients if gravity 
is Lorentz invariant (see (2.68)). 

And generalizing the geodetic precession (2.32), a gyroscope in orbit will undergo the 
secular precession; 

dS/dt = Va(c, + Vfec 7 ) ^ x S 2.67d 

R 3 


along with a non-secular variation in spin axis in non-metric theories of gravity; 
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2.67e 



s vv 


The placement of an accurate gyroscope in low orbit around the Earth would permit observing 
both the precession (2.66) and possible anomalous contributions (2.67a-e). <7) 

Summarizing these results, there are the following constraints from various observables; 


Equation 

Constraint 

2.6 

c, + c 2 + c 3 + 1 = 0 

2.9 

2(c 7 + c, + c 9 + c, 0 ) - (c 4 + c 5 + c 6 ) - 3/2 = 0 

2.14 

c 7 + c, + c, + c l0 = 0 

2.16 

c 7 + c„ = 0 

2.19 

2(C| ■ Cj) + c 7 + c 9 - Cg - C|q — 0 

2.23 

2c, + Cj + 2 ( 0 , + c 9 ) = 0 

2.25 

c 7 + c, - 1 = 0 

2.28 

2c 4 + c 3 = 0 

2.51 

+ c 5 + c 8 + 3/2 = 0 

2.58a 

1 + Cj + 2cj + (c s + 2c s )/3 = 0 

2.58b 

9 t "*■ c,o - (Cj + 2c 6 )/3 = 0 


All of the above observables are related to preferred frame (w dependent) effects. Two additional 
observables are independent of preferred frame effects; 

2-33 c, + Vic-, - 3 

2 42 c 7 * 1 - a - 2r 

The first eleven constraints are redundant, but self consistent; they are all fulfilled by the 
relationships; 
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c 4 = 0 


(gauge or coordinate choice) 


2.68a 


c 6 = -3/2 + c 4 

2.68b 

c, = -2 c 4 

2.68c 

c 3 = 1 + c, 

2.68d 

c 7 = -2 - 2c, 

2.68e 

c, arbitrary 

2.68f 

C « = * C 7 

2-68g 

c 9 = 1 - c 7 

2.68h 

c io - c 7 - 1 

2.68i 

c 7 arbitrary 

2.68j 


These constraints on the c, guarantee that no physical effects proportional to w or (w) 1 will be 
observed in a system moving with respect to the cosmos, when that system is observed from its 
own rest frame. 

Geodetic precession, as observed in the lunar orbit, then provides a constraint between 
c, and c 7 (2.33)/ ,0) And absence of anomalous changes in the radius or period of satellite orbits 
as they follow the eccentric orbits of a planet around the Sun, specified c 7 in terms of the well 
measured speed of light function in a gravitational environment, and the gravitational time 
dilation of non-gravitational clocks (2.42). For a = 1, T = -2 (2.33) and (2.42) yield the specific 
coefficient values; 
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relativity theory (in an appropriate gauge). 


In particular, the coefficients which determine the observables related to the dragging of 
inertial frames — (2.62), (2.66) and (2.67) — are then numerically uniquely specified; 

c a + c * = *2 -2c, - c 7 = -8 2.70a 

c j - c g + 2cj/3 = 0 2.70b 

From (2.70a) it is seen that the normal frame dragging coefficient (c, + c g ) is predicted from the 

geodetic precession coefficient (2.33) if in addition all preferred frame effects are observed to 
be absent. 

III. Constraints on the Non-Linear Terms in the Gravitational Equation of Motion. 

The last term of (2. 1 ) represents the intrinsically non-linear, post-Newtonian corrections 
to linearized gravity, as well as acceleration dependent terms from linearized gravity in which 
the acceleration of various bodies have been set equal to their Newtonian values, thereby 
producing non-linear contributions to the equation of motion. 

Expressing the non-linear part of (2.1) as a phenomenological expansion continues the 
philosophy of this investigation; 




jk 


jk 


3 

*ik 


c Vc r ‘ k 

"T 12 T 

« r jk 




M3 


+ C 


vV«j 


14 


+ c 


r mm 

i5 — 5 — r ik 


V«J _ f Jk f Jk 


‘Ik 


L r u 


> J 


3.1 


33 


The sum over j and k can include body i, itself, in some cases. The only restriction is that the 
two body indices on one interbody vector can not be the same. A strong constraint on the type 
of terms that can be included in the phenomenological expansion (3.1) results from the observed 
isotropy of the Newtonian interaction. In (3.1) one can let body j or k become a distant body 
(spectator body) which then effectively renormalizes the strength of the two body interaction. 
(3.1) becomes in this spectator limit; 




^ c r m A 

^ 15 — 7“ 


j 


M) 


RR - R R 
* • 8 8 


R 


3.2 


If the c 15 term was not the exact difference between two terms which canceled in the spectator 
limit, it would yield a Newtonian interaction between i and j which was dependent on the angle 
between r fJ and the spectator location R,. Considering the galaxy as the spectator body to the 
solar system (M^R, - 10" 6 ), the observed isotropy of the Newtonian interaction to a part in 10' 3 
strongly constrains (at the 10' 7 level) the type of terms which can appear in (3.1). (,,) 

We first consider observables within the two body problem but which include non-linear 
contributions from (3.1). In this case (3.1) simplifies to; 


■ c. 


m i% 


+ ( C 12 ~ C I3 C u) 


mjm.r.j 


3.3 
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The solar system mass is so dominated by the Sun’s mass, the second term in (3.3) plays no role 
in any Sun-planet observables; c u is the only non-linear coefficient measurable to good accuracy 
from observations of post-Newtonian corrections to planetary orbits about the Sun. From ( 3 . 3 ) 
and the rest of the equation of motion ( 2 . 1 ) planets or test bodies orbiting the Sun are governed 
by the equation of motion (c 4 = 0 gauge); 


a = 


Mr r 


1 + c,v 2 


Mr ■ v v M 7 ? 

+ c. + c,, 


3.4 


(Observable 15) For a circular orbit (3.4) yields a modified Kepler’s third law relationship 
between orbital period and orbital coordinate radius; 



1 + < c . " c „> 


M 


3.5 


Observations by radar ranging between the inner planets will give the weak constraint; 

c i - c M = -3 3 g 

(Observable 16) For eccentric orbits (3.4) is most transparently solved in the u(0) » l/r(0) 
representation, in which the orbital solutions are given by; 

u( 0 ) = u 0 ( 1 + e cosa 0 } 3 -j 

with 

a * 1 - (c, + c 7 - Viculmu* 3 g 
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(3.7) represents a precessing elliptical orbit with post-Newtonian contribution to orbital precession 
of (per orbital period); 

50 * ( 2 c, + 2c 7 - c,,) 7 tmu 0 3 9 

(3.9) can be best constrained by observation of Mercury’s perihelion precession;* 11 * 

2c, + 2 c 7 - c,, = 6 3.10 

Taking the difference between (3.10) and (3.6) yields a relationship between the two key 
motional coefficients; 

c, + 2c 7 = 9 3 u 

(Observable 17) We now turn to observables which intrinsically require three separate mass 
elements in their construction. The gravitational (passive) to inertial mass ratio of a celestial 
body (particularly the Earth) is an observable which can be measured to good precision. (,J U ,5) 

Considering a celestial body as a gas of particles in internal equilibrium, the gravitational 
self energy contributions to this ratio (gravitational mass/inertial mass) are calculable from ( 2 . 1 ) 
plus the non-linear interaction given by (3.1). Let a body of mass elements m, be accelerated by 
external bodies M k which produce a Newtonian acceleration field; 

* " £ 3.12 

k 

The weighted sum 

m, a, , a, = ij( internal) + a 
1 

is performed over the mass elements of the celestial body, and the tensor and scalar virial 
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£ m.v, 2 - ft £ 
' u 



= 0 


The body’s collective acceleration is then given by; 



3.13b 


3.14 


The constraint required for the gravitational to inertial mass ratio of Earth to be one as observed 
in lunar laser ranging is then; (l6l7) 

3c, - c 7 - 5c, , - c, 2 - 6c, 3 - 2c u - 4c, 3 = 0 3.15a 

Additionally there are the two weaker constraints required for that ratio to be one for a rotating 
celestial body; 

-c? + c„ - c, 2 - 2c, 4 - 4c, , = 0 3.15b 

and 
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3.15c 


C,5 = 0 

(Observable 18) The active gravitational mass of a celestial body can also be evaluated from the 
phenomenological expansions (2.1) plus (3.1). For a gas model of a celestial body, the active 
mass is given by; 
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3.16 


in which R is the unit vector toward the test body measuring the active mass. If the active mass 
is to equal the energy content (first line of (3.16)) of the body, three constraints are required; 
c j + Ch - 2c 12 + Vi = 0 3.17a 

C|4 + 2c 13 - C| 0 - c l3 = 0 3.17b 

and 


c « ■ 3c u = 0 3 17 C 

Unfortunately, there are no solar system observations accurate enough to enforce these 
constraints. 

(Observable 19) The perturbations of the orbit of an earth satellite are of interest. The Earth’s 
orbit about the Sun and the satellite’s orbit about the Earth are assumed to be nominally circular 
and coplanar. The relevant perturbations of the satellite relative to Earth are then; 
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5a = -(2c, + c 2 ) 


mru -v , . mr • vu 
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The resulting radial perturbations of the satellite orbit are then given by; <, *• 1 ’ , 


3.18 


5r(t) = 


cos(g) - Cl)t 


(1 - n/co)(co„ " (® - Cl) 2 ) 


-(2c, + c 2 + 2c 7 + 2c 9 ) JH.uv 

r 2 

(3c, 3 + 2c,, + c, 4 + c, 3 - c, + c 2 + 2c 7 ) 

Cl , 3 . 

"® (C » +C ”" C * + 4 ^ ) 


m M 
T R 2 


3.19 


R and r are the positions of Earth from the Sun and satellite from the Earth, respectively, v * 
Rfl and u = ra> are the speeds of the Earth relative to the Sun and the satellite relative to the 
Earth, respectively. Cl and to are the orbital angular frequencies of the Earth and satellite, 
respectively. (0 o is the frequency of perigee of the satellite orbit. M and m are the masses of 
the Sun and Earth, respectively. 

The two largest contributions to (3.19) are dimensionally larger than the experimental 
accuracy with which lunar laser ranging (and ranging to low earth satellite orbits) data is fitted 
by the general relativistic model; hence there are the two constraints; 

2c, + (h + 2c 7 + 2c, = 0 3.20a 

and 

3c, j + 2c„ + c I4 + 5c ts /4 - c, + Cj + 2c 7 = 0 3.20b 
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The first of the above constraints is the same as (2.23). 

All of the constraints presented in this section on the non-linear gravitational interaction 
are satisfied with the coefficients in (3.1) taking the general relativistic values; 


c„ = 4 

3.21a 

C I2 = 1 

3.21b 

c„ = -7/2 

3.21c 

c„ « -1/2 

3.21d 

C|3 = 0 

3.21e 


along with the first ten c, taking the values given in section n, which are also the general 
relativistic values. However several of the constraints of this section are too weakly specified 
by observation, so that all five of the non-linear coefficients are not determined by solar system 
observations alone. We do obtain the weak constraint (3.11) between c, and c 7 though through 
consideration of these observations involving non-linear aspects of the gravitational interaction. 
IV. Experimental Accuracy of Constraints on Observables 

The models used to fit solar system observations are not as general as (2.1); they will 
generally assume at least that gravity is a metric field based interaction. Therefore, a number of 
the observables calculated in this study will not explicitly have been fitted in past analysis of 
observational data. So in this section we will sometimes have to infer what experimental 
accuracy would be if existing data were fit for the observables within our model. 

We give an overview of the experimental accuracy of the various observable constraints 
derived in this study. We first consider two very high accuracy observations which give us the 
strong constraints; 
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c 7 + c g + c 9 - Cg = o ± 10' 7 


4.1 


and 

2(c, - Cj) + c-j + c 9 - Cg - C|Q = 0 i 10^ 4.2 

(4.1) results from the observed isotropy of the Newtonian gravitational interaction between the 
mass elements of the Sun: if (4.1) was not fulfilled, the Sun’s spin axis would have processed 
out of solar system alignment over the past 4.5 10 9 years. 00 (4.2) is the result of a possible self- 
acceleration of a spinning celestial body; 


a 



+ C 7 + c 9 - c , - C, 0 ] 


4.3 


U is the body’s gravitational self energy, Cl is its rotational angular frequency. Application of 
this to the Earth, as observed with lunar laser ranging data, leads to the constraint (4.2). < * 0> 

The constraint (2.23) is fairly strong; 

2c, + Cj + 2c 7 + 2c 9 = 0 ± 10‘ 3 44 

and results from absence of anomalous range oscillations in low Earth satellite orbits (3.19) as 
measured by laser ranging. 

The absence of gravimeter anomalies on Earth larger than about 10‘ 7 gal give the 
moderately strong constraints;* 20,20 


C4 + Cj + c tf + 3/2 = 0 ± 10' 3 

4.5a 

1 + 2c 2 + 2c, + (c, + 2c s )/3 = 0 ± 10' 2 

4.5b 

c, + c 10 - (c, + 2c a )/3 = 0 ± 10' 2 

4.5c 


The constraint (2.16) is quantified by the upper limits on anomalous secular or 
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semisecular changes in Earth’s orbital period; 

c 7 + c, = 0 ± 10' 3 46 

While (2.25) is required to suppress longitudinal anomalies in the inner planet position greater 
than a few tenths of a kilometer, and as measured by radar ranging; 

c 7 + c, = 1 ± 10- 2 4 ? 

Finally (2.6) is quantified by the absence of anomalous period changes in earth satellite orbits; 
c i + c 2 + c 3 + 1= 0± 10' 3 4 g 

The solution (2.68) for the motional coefficients is therefore accurate to about 1 percent, with 
some of the relationships much stronger. 

The geodetic precession observation of the Moon’s orbit is, however, only about 10 
percent accurate. 00 ' So without the use of other observables involving non-linear aspects of the 
gravitational interaction, the dragging of inertial frames proportional to c, + c, can be predicted 
to only about 10 percent accuracy. In fact, the inclusion of the non-linear observables in section 
m does not improve the accuracy with which c, and c 7 can be determined. This is because a 
number of new coefficients must be introduced to express the general non-linear interaction (3.1); 
and the number of independent observables of high accuracy in the solar system does not 
increase sufficiently to both constrain the new coefficients in (3.1) as well as c, and c 7 . 

V. Conditions on the Equation of Motion for Momentum Conservation. 

The phenomenological gravitational many body equation of motion (2.1) as it stands in 
its generality does not possess a conservation law for total momentum of an isolated system of 
bodies. However, certain constraint relationships among the coefficients c, will indicate a 
momentum conservation law. Let an isolated system’s total conserved momentum be given by 
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an expression of the general form; 


E fo + p,"W + p,E — + p,E W 

1 J r u j r„ 


5.1 


with p ) I 3 being three dimensionless constants. Then dP/dt = 0 yields; 
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Using (2.1) and (3.1) to express the first term in (5.2), and elsewhere using the Newtonian 
acceleration field to represent a, in post-Newtonian terms of (5.2), i.e.; 

a, * E m i 

k 

some of the coefficients c, in (2.1) and (3.1) are then constrained to relationships involving the 
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three constants p,. 


Ci - c 3 = *Pi + p 3 
c 4 - c s = -3p 3 

C 7 ” C IO = 2p t + p 2 - Pj 

c« - c, = p 2 + p 3 

C ll ' C 12 + C I3 = P2 

C|4 + 2c i3 = p 3 

If a center of energy exists which is to move most generally at constant velocity; 
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5.3a 

5.3b 

5.3c 

5.3d 

5.3e 

5.3f 


5.4 


with 


dR/dt - P/E 


and therefore 


d J R/dt J = 0 


then 


Pi = -P 2 - -Pj = 1/2 


5.5 
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and the constraints (5.3) are then even more restrictive. Not surprisingly, the momentum 
conservations constraints (5.3) automatically lead to the vanishing of the observable (2.19) 
associated with the possible self-acceleration of the center of mass of a gravitational system. 
VI. Conditions on the Gravitational Equations of Motion to be Metric Field Derivable. 

The phenomenological gravitational many body equations of motion (2.1) plus (3.1) are 
more general than what can be obtained from a metric theory of gravity. In metric theories there 
are several potential functions — a symmetric g^r ,t) — which for test bodies yield the equation 
of motion; 

d(9L/9v)/dt - 9L/9F = 0 6.1a 

with 


L(r,t) = - y/ giiv (F,t)dx7dtdx7dt 


6.1b 


A general expansion for the metric field potential contains seven free parameters; 
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6.2b 
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8,j - -(1 -M,U)5„ 


6.2c 


with 


U(7,t) = y , 

J 



6.2d 


Four coordinate gauge choices were available to set some of the metric terms generally equal to 
zero: the spatial g,, is kept diagonal and isotropic by a spatial coordinate transformation; 

7 + m^-r^/lr-rj 6.3a 


while an acceleration dependent term in g^ is generally eliminated by a time transformation; 

* t+ xE m k (r-r k )*v k /l 7_7 kl 6.3b 

k 


Since there are fewer free parameters in (6.2) than there are coefficients c, in (2.1) plus 

(3.1) , constraints on the c, are necessary in order that the equation of motion is derivable from 
a metric field of general form. Using (6.2) in the Euler-Lagrange equation (6.1) yields an 
equation of motion whose terms can be compared to die phenomenological equation of motion 

(2.1) plus (3.1). The constraints on the c, necessary for a metric field based equation of motion 
are then; 

C| = -W2 6.4a 
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Cj = -(M, + M 4 ) 

6.4b 

c 3 = M 4 - M,/2 

6.4c 

II 

O 

in 

II 

o 

6.4d 

c* = -30V2 + M 4 ) 

6.4e 

c 7 = 2 - Mj 

6.4f 

c, = -(M 3 + M 4 ) 

6.4g 

c, — M, - 1 

6.4h 

C|o = M 3 - M 4 - M2 

6.4i 

c„ = Mg - M 3 + 1 

6.4j 

c 12 = Mj/2 
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Several of the observables are determined by the fulfillment of the above metric constraints. For 
example, (2.25) and (2.28) are automatically fulfilled if (6.4) hold. Also the anomalous 
gyroscope precession (2.67a) vanishes if gravity is derivable from a metric field. 

VII. Invariance of Motional Observables Under Coordinate (Gauge) Transformations. 

The N body coordinates r, used in (2.1) can be combined to form N new, independent 
coordinates associated with the motion of the N bodies. The equations of motion (2.1) and (3.1) 
would of course be altered under these coordinate transformations. In particular, if the 
transformation; 
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7.1 


r . * 7 ( + £E "\ r A 

fari 


is made, the ten c, in (2.1) are altered by; 


c , C, + 5 

7.2a 

Cj Cj - 2 % 

7.2b 

c 3 -> c 3 + % 

7.2c 

c< — > c 4 - 31; 

7. 2d 

c 5 -* c 3 + 6^ 

7.2e 

c # c 6 - 3^ 

1.21 

c, -* c 7 - 2 % 

7-2g 

c, -> c, + 2 % 

7.2h 

c, c, + 2^ 

7.2i 

C 10 C IO - 2^ 

7.2j 


The form of changes also. The coordinate speed of light function is altered, a time derivative 
of (7.1) giving in first order approximation; 

One is therefore always free to choose a gauge or set of coordinates which sets one of 
= However, the observables which were derived in section II are all invariant under this 
transformation. Setting c 4 = 0 is the conventional choice of coordinates, as this leads to a simple 
scaling relationship between the spatial global coordinates and the local proper spatial coordinates 


c(r,t,6) 


c(r,t,e) + ^|U - £ m k ((r - r k ) • / 
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measured, in particular, by light ray propagation; 

5r p = (1 + U(r,t))5r 7 3 

It is important to note that the speed of light function is isotropic (independent of c P = 0) 

c(F,t> = 1 -i£ ny|7-7| *r£ m.«7-7.) •«>’/ |7-7| J 

J J ] 

in the same coordinate system that ( 7 . 3 ) is valid. 

VIII. Gyroscope Precession. 

In our general, phenomenological model for the gravitational equation of motion of 
matter, we derived several possible contributions to the precession rate of the spin axis of an 

orbiting gyroscope. There are the two well known secular terms which are non-zero in general 
relativity; 

d5/dt = ^(c, + Vic-,) m (r x dr/dtVr 3 x 5 g 1 


+ ( C j + c^/8 (J - 3J ?f)/r 3 x S 


( 8 . 1 ) is the geodetic precession contribution, while ( 8 . 2 ) is the "frame dragging" precession due 
to the Earth’s spin angular momentum 7. (In general relativity c, = 1 , c 7 = -c, = -c g = 4 ). 
There is an additional secular precession in non-metric theories of gravity; 


ds/dt = - ± L - C# + l c ,) 

16 1 3 J a 0 3 (l-e 2 ) 3 ' 2 


8.3 
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a »’ e P b®* 11 ? th* gyroscope orbit’s semi-major axis, eccentricity and unit polar vector, 
respectively. 

Non-secular precessions of sufficient magnitude to be measured are also of interest. In 
gravitational theories which have a preferred inertial frame (absence of Lorentz invariance), there 
are some possible changes in the gyroscope spin axis; 

«S =Va(4 - 3^ - c, - 2c, + 2c, + c 2 + 2c 4 + c 5 ) V • wS 

+ (2 + 2c, + c 2 + c 7 + Cg)/4 (w x V) x i g 4 a 

+ (c 7 + c, - 2 - 2c, - c 2 )/4 (S ■ Vw + S wV) 


and the additional precession rate (not integrable in closed form); 

dS/dt = %(2c 4 + Cj) MS *wR/R 3 g 45 

w is the velocity of the solar system with respect to the cosmological preferred inertial frame, 
and V = dR/dt is the velocity of the gyroscope with respect to the Sun, which has mass M. 
These precessions (8.4a,b) are absent in general relativity, but dimensionally are of magnitude 
80 mill iarc seconds and have an annual period. 

Finally there is a small non-secular variation in gyroscope spin axis which only exists in 
non-metric theories of gravity; 

S = ^(Vic, - c, - 2cJ3 - 1) S • V V g 5 
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which is dimensionally only of order 2 milliarcseconds. 

It should be pointed out that the geodetic precession (8.1), though dominated by its earth 
orbiting contribution (-3 10 3 mas/yr), also has a contribution (-13 mas/yr) from the orbital 
motion of the gyroscope about the Sun. 

As was shown in previous sections, all the coefficients in the motional, post-Newtonian 
gravitational interaction — including the gravitomagnetic coefficients — are measured to some 
precision by other observations in the solar system which do not involve gyroscope precession. 
Therefore the precession rates can be predicted with accuracies of about 10 percent, or perhaps 
slightly better. One could then view the gyroscope precession observations as a way of 
substantially improving the measurement of the coefficient combination in (7.1) — c, + Vic v 

If the gyroscope precessions are found to be different than predicted, we believe this 
would present a major crisis for gravitational theory. The model (2.1) is quite general, especially 
with regard to the motional, two-body post-Newtonian interaction which includes 
gravitomagnetism. Given two mass elements located at relative coordinates =7, -7j, and 
each having velocity v„Vj, respectively, we have found it impossible to generalize a post- 
Newtonian acceleration expression beyond the form given by (2.1). Yet this would seem to be 
necessary if redundant observations could not all be predicted by the nine coefficient model (2.1). 

Since the spin angular momentum of the gyroscopes to be placed into orbit is motional 
angular momentum — not quantum spin angular momentum — there would appear to be no way 
to explain unpredicted precession by means of some anomalous coupling of gravity to spin. For 
macroscopic gyroscopes, spin angular momentum is simply a manifestation of particle motion, 
superimposed over all the matter in the gyroscope. 
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